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A  STUDY  OF  ADAPTIVE  DETECTION  OF  RANGE-DISTRIBUTED  TARGETS 


I.  INTRODUCTION 

A  recent  paper  by  Bose  and  Steinliardt  [1]  significantly  advanced  the  theory  of  space-time  adaptive 
detection  of  a  desired  signal  in  the  presence  of  correlated  Gaussian  noise.  Space-time  adaptive  detection 
is  concerned  with  detecting  a  signal  of  partially  known  direction  and  doppler  and  unknown  strengths 
across  time  (or  range)  in  the  presence  of  Gaussian  noise  whose  covariance  matrix  is  unknown.  Their  work 
generalized  much  of  the  adaptive  detection  research  [2-6]  of  the  last  twenty  years.  A  variety  of  application 
problems  can  be  solved  using  their  formalism.  They  did  this  by  generalizing  the  desired  signal  waveform 
to  be  a  rank  one  waveform  (in  matrix  form)  that  is  uncertain;  i.e.  only  partially  specified.  More  precisely, 
they  modeled  the  space-time  steering  vector  as  a  Kronecker  product  of  two  vectors,  each  of  which  is 
unknown  but  is  known  to  lie  in  a  known  subspace. 

One  limitation  of  their  approach  is  that  the  dimensions  of  these  subspaces  are  constrained  with  respect 
to  the  dimensions  of  the  space-time  array  and  the  number  of  snapshots  (data  vectors)  associated  with  the 
space-time  array.  These  constraints  were  imposed  because  in  the  derivation  of  the  Generalized  Likelihood 
Ratio  Test  (GLRT)  detector  it  is  possible  that  a  certain  estimated  covariance  matrix  could  become  singular 
which  would  significantly  complicate  the  derivation.  Hence  constraints  were  imposed  to  eliminate  the 
possibility  of  generating  singular  matrices.  We  introduce  what  we  term  a  Modified  Generalized 
Likelihood  Ratio  Test  (MGLRT)  detector  in  order  to  generate  a  useful  detection  structure  for  the  case 
when  the  estimated  covariance  matrix  is  singular. 

In  this  paper,  we  solve  a  problem  of  practical  importance  which  cannot  be  solved  by  the  formalism 
of  [1].  Consider  the  problem  of  adaptively  detecting  a  range  distributed  target  or  targets.  Suppose  the 
desired  target  or  targets  could  be  spatially  distributed  across  the  entire  range  interval  whose  returns  are 
to  be  used  as  the  input  data  block  (sensor  x  range)  from  which  the  adaptive  detector  is  derived.  Assume 
that  the  scatterers  remain  in  their  respective  range  cells  during  a  coherent  processing  interval  (CPI);  i.e. 
the  scatterers  do  not  range  walk  during  a  CPI.  The  scatterers  returns  may  be  uniformly  doppler  shifted 
due  to  target(s)  motion.  If  one  uses  a  single  point  target  adaptive  detector  (such  as  Kelly’s  [2]  or  the 
Adaptive  Match  Filter  Detector  [4,6])  to  detect  individual  point  scatterers  across  the  range  interval  of  the 
data  block,  significant  detection  performance  degradation  could  occur  in  each  range  cell  where  a  target 
scatterer  is  present  due  to  the  desired  signal  contamination  of  the  secondary  data  (data  where  it  is 
incorrectly  assumed  a  desired  signal  is  not  present)  [7].  Hence  the  total  adaptive  detection  performance 
(detecting  the  presence  of  any  target  or  targets  across  the  range  interval  of  interest)  could  be  significantly 
degraded. 

We  derive  the  MGLRT  associated  with  adaptively  detecting  a  range  distributed  target  or  targets.  It 
is  pointed  out  that  our  methodology  is  not  restricted  to  detecting  a  single  target  that  is  range  distributed 
(resulting  from  using  a  high-range  resolution  (HRR)  radar  waveform)  but  could  also  be  used  to  detect 
a  formation  of  targets  with  the  same  velocity  that  are  spatially  distributed  in  range  (say  from  a  low  range 
resolution  radar  waveform).  In  fact  the  problem  of  signal  cancellation  that  occurs  when  one  performs 
single  point  target  adaptive  detection  can  be  significantly  abated  by  using  the  proposed  methodology: 
instead  of  trying  to  detect  individual  targets  and  suffering  the  potentially  large  detection  losses  due  to 
signal  contamination,  one  should  try  to  detect  the  presence  of  the  target  set  within  the  input  data  block. 
A  benefit  of  doing  this  is  that  there  will  be  an  effective  integration  gain  of  target  returns  which  enhances 
the  target  set’s  detectability.  Hence  what  was  thought  to  be  a  source  of  performance  degradation  (signal 
contamination)  can  actually  be  a  source  of  performance  enhancement.  Obviously,  the  range  resolution 
of  the  target  could  suffer  via  this  methodology  but  at  least  a  detection  may  not  be  lost. 
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The  outline  of  this  paper  is  as  follows.  The  MGLRT  for  a  range  distributed  target(s)  is  defined  and 
derived  in  Section  II.  Expressions  for  false  alarm  and  detection  probabilities  are  derived  in  Section  III. 
Unfortunately,  the  MGLRT  for  a  single  data  block  is  not  a  constant  false  alarm  rate  (GEAR)  detector, 
i.e.  the  test  statistic  is  not  independent  of  the  external  noise  environment.  In  Section  IV,  it  is  shown  with 
the  use  of  multiple  data  blocks  how  to  construct  a  GEAR  detector.  In  the  preceding  development,  it  has 
been  assumed  that  the  desired  steering  vector  was  known.  In  Section  V,  we  derive  the  MGLRT  for  range 
distributed  targets  when  it  is  assumed  that  the  desired  steering  vector  is  unknown  but  lies  on  a  known 
subspace.  Some  simulation  results  are  presented  in  Section  VI.  In  Section  VII,  the  single  point  target 
MGLRT  with  assumed  contamination  of  the  secondary  data  is  derived.  However  the  resultant  MGLRT 
is  useless  since  it  is  shown  that  none  of  the  secondary  data  is  used.  This  was  an  attempt  to  solve  the  long 
standing  problem  of  how  to  detect  a  point  target  in  contaminated  secondary  data.  The  MGLRT  formalism 
turned  out  not  to  be  an  effective  solution. 

11.  MGLRT  FOR  RANGE  DISTRIBUTED  TARGETS 

A.  Preliminaries 

Data  is  collected  from  N  sensors  for  a  radar  that  has  a  transmitted  waveform  with  a  fixed  pulse 
repetition  interval  (PRI).  For  example,  the  V  sensor  inputs  could  consist  of  an  array  of  antennas  and  time 
taps  (with  the  time  delay  equal  to  the  PRI)  of  these  antennas.  A  data  vector  consisting  of  K  elements  is 
collected  from  each  sensor.  The  data  associated  with  the  Ar’th  element  of  each  sensor  data  vector  is 
assumed  to  represent  data  coming  from  the  A:’th  range  cell.  Assume  that  the  target  (or  targets)  is  spatially 
distributed  across  all  of  the  K  contiguous  range  cells  in  some  fashion  (some  range  cells  could  have  no 
target  scatterers).  The  returns  for  the  K  contiguous  cells  on  each  of  the  N  sensors  could  consist  of  the 
target  returns  plus  interference  or  interference  by  itself.  The  interference,  assumed  to  be  strict-sense 
stationary,  could  come  from  a  variety  of  sources:  e.g.  clutter,  jamming,  and  system  noise.  Inputs  are 
complex  valued. 

It  is  pointed  out  that  if  we  had  assumed  that  only  a  specified  m  contiguous  range  cells  out  of  the  K 
range  cells  could  contain  target  returns  and  K  —  m  >  N,  then  the  methodology  of  [1]  could  be  applied. 
However,  for  our  development,  it  is  reasonable  to  assume  that  we  do  not  know  which  of  the  K  range  cells 
contain  the  target  or  targets.  Thus  we  assume  that  the  target  returns  could  occur  anywhere  in  the  K  range 
cells. 

We  desire  to  construct  a  hypothesis  test  which  distinguishes  between  the  signal  plus  interference 
hypothesis  (Hi)  and  the  interference  only  hypothesis  (Hq).  Gonsider  the  two  hypotheses: 

Hi  :  +  a/s  ,  k  =  1,  2,  . . K 

where  x*,  z^.,  s  are  the  interference  vector,  received  vector,  and  desired  signal  steering  vector, 
respectively,  all  of  length  N  associated  with  the  ^-th  range  cell  of  the  N  sensors.  The  desired  steering 
vector,  s,  is  normalized  so  that  s“s  =  1  where  H  denotes  conjugate  transpose.  For  examples,  the  steering 
vector  for  the  space/time  adaptive  processor  (STAP)  is  the  Kronecker  product  of  desired  angular 
(direction  of  arrival)  steering  vector  and  the  desired  doppler  shift  steering  vector.  The  complex  scalar 

where  *  denotes  complex  conjugation  is  the  unknown  complex  signal  amplitude  of  the  desired  target’s 
scatterer  in  the  ^’th  range  cell  (we  conjugated  for  notational  purposes).  Eq.  (1)  can  be  rewritten  more 
succinctly  as 
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(2) 


//„:  Z  =  X 
//,  :  Z  =  X  +  sa» , 

where  Z  and  X  are  denoted  as  the  N  X  K  input  data  matrix  and  interference  matrix,  respectively,  a  = 
(<2i,  (22,  and  T  denotes  transpose.  The  input  interference  vectors,  X/.  (k  =  1,  2,  ^  are 

assumed  to  be  independent  and  identically  distributed  (i.i.d.)  zero-mean  complex  circular  Gaussian 

vectors.  The  covariance  matrix  of  each  vector  is  given  by  an  N  x  //complex  Hermitian  matrix,  R, 

which  is  assumed  unknown.  Assume  /iT  >  A^.  For  AT  >  A^and  R  positive  definite,  the  rank  of  Z  equals 
A/' with  probability  1  [10]. 

It  can  be  shown  (see  e.g.  [8])  that  the  probability  density  function  (pdf)  for  Z  under  each  hypothesis 
can  be  written  as 


:  p^{Z-  R)  =  -S—  exp[-  rr{/?*ZZ"}]  5(-f>”Z) 

//,  :  p,{Z-  R,  a)  =  exp[-  Tr{R\Z  -  sa")(Z  -  sa«)"}]  5(^2  (Z  -  sa")). 


(3) 


where  c  is  the  pdf  normalization  constant,  |A||p  denotes  the  positive  determinant  of  R  (the  product  of 

the  positive  eigenvalues  of  R),  Tr{  }  is  the  trace  of  the  argument,  R^  denotes  the  Moore-Penrose 

generalized  inverse  of  R,  6  is  the  product  of  Dirac  delta  functions  of  the  matrix  elements  of  the 
argument,  and  4>2  will  be  defined  shortly.  The  generalized  inverse  is  defined  as  follows.  We  can  write 

R  =  ^>AT>"  where  is  the  N  x  N  matrix  of  eigenvectors  of  R  and  A  is  the  diagonal  matrix  of 
eigenvalues  of  R.  Write  A  =  diag  (X,,  . . .,  0,  . . .,  0}  where  X„,  >  0,  m  =  1,  2,  . . .,  Mand 

M  <  A/and  A-"  =  diag  {X,'',  X2',  . . .,  X^^',  0, . . .,  0} .  Then  A*  =  $A^4>"  [9],  Set  $  =  («i>,  ^>2)"  where 
the  {N  —  M)  X  N  matrix,  $2,  is  associated  with  the  null  space  of  the  columns  of  R. 

In  the  adaptive  detection  literature,  it  is  almost  always  assumed  that  R  is  nonsingular  so  that  R*  =  A  ' 

11^ lip  =  II ^11  where  H-l  denotes  the  determinant  of  .R.  However,  we  will  find  in  our  formulation 
that  this  assumption  need  not  hold.  We  will  find  that  when  we  formulate  a  modified  maximum  likelihood 

(MML)  solution,  for  R  under  that  /?mml  be  made  singular  for  a  certain  subspace  of  a. 
The  GLRT  is  formulated  as  follows. 


sup  p^(Z;  R,  a)  //, 

^ -  >  t.  (4) 

sup  P(,(Z;  R)  < 

R  ^0 

where  R  is  the  estimate  of  R  under  each  hypothesis.  It  will  be  found  for  our  problem  (where  po  and  /?, 
are  defined  by  (3))  that  sup  p^{Z\  R,  a)  is  unbounded.  The  pdf  p^  is  infinity  as  the  range  space  associated 

/?,a 

with  the  columns  of  R  overlaps  the  null  space  associated  with  the  columns  of  Z  -  sa".  The  null  space 
of  Z  -  sa”  is  non-empty  if  we  can  choose  an  a  such  that  Z  -  sa“  has  rank  less  than  N.  If  Z  has  full 
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rank,  we  can  always  find  an  a  that  drops  the  rank  of  Z  —  sa"  to  -  1.  For  example,  set  a  =  Z”s. 
Equivalently,  p^  approaches  infinity  as  {R,  a)  approaches  (with  respect  to  some  distance  measure)  a 
member  of  the  following  set,  S 

S  =  {R,  a  \  p{R)  =  N  -  I,  R  =  (Z  -  sa“)  (Z  -  sa")"} 

where  p(*)  denotes  the  rank  of  the  argument.  Here,  we  see  S  consists  of  the  subspace  of  a  where  (Z  — 
sa")  (Z  -  sa")"  is  a  singular  matrix.  Thus  it  would  seem  that  a  useful  GLRT  for  this  problem  does  not 
exist  since  sup  p^  is  unbounded.  As  Grenander  [13]  points  out,  in  certain  ML  problems  the  parameter 

/?,a 

space  is  too  large.  He  proposes  salvaging  the  method  by  performing  the  maximization  over  a  non-trivial, 
constrained  set.  This  method  is  called  the  “Method  of  Sieves.”  Using  the  Method  of  Sieves  [13],  we 
can  restrict  (R,  a)  to  a  reasonably  subspace  such  that  a  maximum  likelihood  solution  exists.  Intuitively 
it  would  seem  that  if  there  were  a  useful  estimate  of  (R,  a),  it  should  come  from  set  S  since  there  exists 
(R,  a)  in  the  neighborhood  of  a  given  member  of  S  which  yield  an  arbitrarily  large  value  for  pi.  Thusly, 
it  will  be  found  that  a  reasonable  estimator  and  detector  can  be  found  for  this  problem  when  we  restrict 
(R,  a)  to  be  a  member  of  S.  However,  on  this  subspace  we  maximize  with  respect  to  R  the  multipliers 

on  the  5  ftmctions  seen  in  Eq.  (3).  Set  m^{Z-,R)  and  /n,(Z;  R,  a)  equal  to  the  multipliers  of  the  6 
functions  for  Po  and  p,,  respectively  seen  in  Eq.  (3).  On  this  subspace  for  cases  of  interest,  sup  m,  (Z; 
R,a)<  oo  since  is  defined  by  (3)  for  all  R  with  rank  greater  than  or  equal  to  one.  It  will  be  seen  in 
the  following  subsection  that  we  will  have  to  further  restrict  the  form  of  R  in  order  to  obtain  a  unique 
ML  solution.  Hence  this  suprema  will  be  used  in  the  likelihood  ratio  test  expression  and  a  finite  test 
statistic  results  (see  [15]  for  other  examples  of  where  the  Method  of  Sieves  is  used  for  estimating 
covariances  matrices).  We  note  for  the  case  when  N  =  I  that  sup  tn^(Z;  R,  a)  is  unbounded.  This  case 

R,a€S 

was  presented  in  [14]  as  a  counterexample  to  the  statement  that  a  maximum  likelihood  estimator  always 
exists. 

When  we  restrict  the  subspace  of  the  unknown  parameters  and  find  the  maximum  of  the  multiplier 
of  the  delta  function  for  each  pdf  of  the  received  data  (under  either  hypothesis)  with  respect  to  the 
unknown  parameters,  we  call  this  the  modified  maximum  likelihood  (MML)  estimate.  When  we  use  the 
ratio  of  the  m,  and  Wq  estimates  we  call  this  the  modified  GLRL  (MGLRT).  As  in  the  case  of  the  GLRT, 
there  is  no  claim  to  optimality  as  a  detector  for  the  MGLRT.  However,  as  in  the  case  of  the  GLRT,  the 
MGLRT  seems  like  a  reasonable  approach  for  finding  a  good  detector.  This  property  of  “goodness”  must 
be  checked  out  via  simulation  and  comparison  with  alternate  approaches.  We  will  see  that  the  MGLRT 
for  this  particular  problem  is  a  useful  detector.  However,  as  we  will  show  in  Section  VII  for  a  different 
detection  problem,  a  useful  solution  does  not  always  result. 

B.  Modified  Maximum  Likelihood  Solution  for  a  Singular  Sampled  Covariance  Matrix 
Let  Z  have  the  following  pdf: 


p  (Z;  R)  =  exp[-  Tr{R*  R}]  6(4>2Z),  (5) 

where  N  >  2,  R  =  ZZ"  can  be  singular  with  at  least  rank  1  and  R  has  at  least  rank  1.  If  p(R)  =  M<N 
then  it  can  be  shown  that  the  ML  solution  for  R  does  not  exist.  However,  we  can  find  a  unique  MML 
estimate  for  R  by  restricting  R  as  follows.  Define  the  set  S^ofN  z  N  matrices: 
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(6) 


=  \R\R  =  $  • 


A  0 
0  0 


where  R  =  is  the  eigenvalue  decomposition  of  R  and  A  is  any  p(^)  x  p(R)  matrix.  Then  we 

will  show,  =  arg  max  mj[Z,  R)',  i.e.,  a  unique  MML  solution  exists  on  Sq. 

Res„ 

We  will  need  the  following  lemma  for  our  development  of  the  MGLRT. 

Lemma  1:  If  we  restrict  R  6  Sq,  the  MML  estimate  solution  for  R  is 


proof:  We  desire  to  solve 


This  is  equivalent  to  solving 


Kul  =  arg  max  m^fZ;  R). 


^MML  =  arg  min  [K  In  ||a||  +  Tr{R*R}]. 

R£S. 


(7) 

(8) 


Rewrite  (8)  as 


=  arg  min  [K  In  \\R\\  +  r/-{R"l>A<l"}] . 


(9) 


This  can  be  shown  to  be  equivalent  to  finding  a  matrix  Rq  such  that 

\  =  arg  min  [K  In  \\Rq\\  +  Tr{Rok}],  (10) 

«oe5, 


where 


R^  =  4>”/?4>, 


(11) 


5i  is  defined  as 


A  0 
0  0 


(12) 


and  A  is  any  p(.R)  x  p(.R)  matrix.  Let  A^  be  the  p(.R)  x  p(R)  diagonal  matrix  of  the  positive  diagonal 
elements  of  A.  Then  (10)  is  equivalent  to  finding 


A^  =  arg  min  In  in  I  +  7rM*A,}].  (13) 

A 


It  is  well-known  that  A^i„  =  \  and  the  lemma  follows. 


□ 
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C-  MGLRT  Derivation 


Using  Lemma  1,  we  find  that 


sup  m^{Z]  R)  =  - £ -  exp  -  rr{(ZZ”)*ZZ"} .  (14) 

||ZZH||; 

We  assume  that  K  >  N  and  the  underlying  covariance  matrix  for  the  A'-length  z  vectors  is  positive 
definite,  so  that  ZZ”  has  rank  with  probability  1,  [10].  Thus 


where  c'  is  a  constant. 
Again  using  Lemma  1 : 


sup  /«o(Z;  R)  = 

R 


c' 

\\zz^V' 


(15) 


sup  ra,(Z;  R,  a)  =  — £ —  exp  -  , 


(16) 


where 


^MML  -  (Z-sa»)(Z-saT.  0^) 

We  note  that  if  R^  has  rank  N  —  I,  then  Tr{R*Rj  =  N  -  I .  Thus 

c" 

sup  m^{Z;  R,  a)  =  - ,  (18) 

"  ||(Z  -  sa«)(Z  -  sa")»||; 

where  c"  is  a  constant.  Thus  the  MGLRT  takes  the  equivalent  form 


min 


\\zz^ _ 

II  (Z  -  saH)(Z  -  saTi, 


fio 


(19) 


Henceforth  for  convenience,  we  shorten  our  notation  for  expressing  the  domain  over  which  the 
minimization  of  the  positive  determinant  is  taken.  It  will  be  assumed  that  the  domain  is  restricted  such 
that  the  matrix  in  the  argument  of  || -Ip  is  singular.  For  example,  we  write  min  I'H  instead  of 

a 

min  II 'Ip.  The  optimal  solution  for  a  is  given  by  the  following  theorem. 
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Theorem  1:  The  solution  for  ~  “  sa")”]]^  is  given  by  a^^,  =  Z"s. 

a 

proof:  Let  be  a  unitary  N  y.  N  matrix  such  that  =  Iq  where  =  (1  0  0  . . .  0)’^.  Let  the  g  - 
decomposition  of  A,Z  =  [L  :  On,k-n]Q  where  L  is  an  x  lower  triangular  matrix  (LTM),  Q  \s  a  K 
X  K  unitary  matrix  and  is  a  N  x  (K  -  N)  matrix  of  zeros.  Now 

Il(z  -  sa«)(Z  -  sa")"!,  ■  U,(Z  -  sa«)!2"e(Z  -  sa")«.4,"|, 

'  IlKi  ;  -  l,a"l[(i  :  0„-,>  - 

where  ao  =  Qa.  Thus  an  equivalent  problem  to  finding  the  optimal  a  is  to  find  the  a^  that  minimizes  the 
far  right  hand  side  (RHS)  of  Eq.  (20).  Equation  (20)  can  be  further  simplified  to  yield 

a^  =  arg  min  ||(L  -  loa”)(L  -  loai")“  +  q;1(,1^||^  (21) 

where  a,  is  an  N-length  vector  of  the  first  N  elements  of  ao  and  a  is  the  sum  of  the  magnitude  squares 

of  the  last  K  -  N  elements  of  ao-  Using  Weyl’s  Monotonicity  Theorem  [9],  the  n’th  ordered  eigenvalue 

of  (L  -  loa“)(^  -  lo^r)”  +  greater  than  or  equal  to  the  n’th  ordered  eigenvalue  of 

(L  -  Iga”)  (L  -  loa,")"-  This  implies  that  (21)  is  minimized  with  respect  to  a  when  a  =  0.  Thus  our 
problem  is  reduced  to  finding  an  AMength  vector  a,  such  that 

a,  =  arg  min  ||(L  -  loa")(L  -  loaJV'llp.  (22) 

a, 

If  L=  (/„„,),  we  will  show  that  a,  =  (/n  0  0  ...  O)"'". 

Set  a,"  =  a^L .  Thus  our  problem  is  equivalent  to  finding  an  N-length  vector  a^  such  that 

a^  =  arg  min  |(4  -  l^a2)LL\I^  -  loa2”)"||/,  (23) 

where  /y^  is  th.t  N  x  N  identity  matrix.  In  order  for  the  matrix  in  the  argument  of  fl  -H^  to  be  singular, 

it  is  required  that  ^  “  Iq®”  singular  which  will  be  if  and  only  if 

32  Ifl  =  1  • 

Using  singular  value  decomposition  (svd),  it  can  be  shown  that 

4  -  V2"  = 

where  Q^  =  [Qj.,,  lol  is  an  N  x  N  unitary  matrix,  is  an  N  x  {N  -  1)  matrix  consisting  of 


(24) 

(25) 
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columns  that  are  orthogonal  to  Iq;  P, 


-La, 


is  an  X  unitary  matrix,  is  an  N  X 


{N  -  1)  matrix  consisting  of  columns  that  orthogonal  to  aj,  and  A,  is  an  A  x  //  diagonal  matrix  with 


diagonal  elements:  X,,  =  Xi2  =  X 


=  X,  =  1 ,  and  X 


0, 


We  can  show 


1(4  -  l„a")LL«(4  -  =  1|A.0»LLHQ,A, 

=  (a^a, )  I  Ji¬ 


lt  is  proved  in  Appendix  A  {Lemma  lA)  that  if  and  <2ii,  are  any  two  N  x  (N  -  1)  matrices  whose 
columns  span  the  space  that  is  orthogonal  to  Ig  then 


(27) 


Hence  we  can  always  find  a  2 , ,  that  is  independent  of  aj  and  the  minimization  of  the  far  RHS  of  (26) 
is  independent  of  the  positive  determinant  term.  Hence,  we  desire  to  find  a^  =  arg  min  a^a^  with  the 

constraint  a  ”4  =  1.  Obviously  a^  =  (1  0  0  ...  0)^ .  Working  backwards,  it  is  straightforward  to  show 
a,  =  L»lo,  io  =  (L:0^^.jv)«4,  and  a  =  V”lo  =  Z”s  =  a^^^-  D 

Collecting  our  results,  we  have  the  following  result. 

Theorem  2:  The  MGLRT  for  detection  of  a  target  (or  targets)  that  are  range  distributed  is  given  by 


|3  = 


|lzz»| _ 

11(4  -  ss»)ZZ“(4  -  ss«)||. 


Ho 


(28) 


We  see  from  the  form  of  /3  that  the  denominator  term  which  was  derived  under  A,  (signal  plus  noise)  will 
be  without  the  additive  desired  signal  vector  because  4  “  ss”  is  the  null  projection  matrix  of  s.  Thus 
(4  —  ss”)Z  will  contain  no  desired  signal. 

We  now  present  a  simple  procedure  for  finding  jS.  A  unitary  matrix  exists  such  that 

=  (0  0  ...  0  l)"^  =  4.  (29) 

Let  ihtQ  -  R  decomposition  of  B^Z  =  (L  :  where  L  is  an  A  x  A  LTM.  It  is  straightforward 

to  show  that  the  test  statistic  /3  defined  by  (28)  reduces  to 
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(30) 


_ 11^^  "II 

1|(/^  -  I„  i:)LL"{I,  -  I,  lo")t 


If  L  =  (/„„),  the  above  reduces  to 


P  -  4. 

Thus  the  procedure  for  finding  the  detection  statistic  is  straightforward: 

1.  Perform  the  Q  -  R  decomposition  of  the  N  x  K  data  matrix,  Z. 

2.  The  NI\r th  element  of  the  N  x  N  triangular  matrix  associated  with  the  Q  -  R  decomposition  is 
the  test  statistic. 

III.  FALSE  ALARM  AND  DETECTION  PROBABILITIES 


A.  Probability  of  False  Alarm 

We  can  derive  the  pdf  of  test  statistic,  (3,  under  the  Hq  hypothesis  by  using  the  equivalent  form 

13  =  _ _  - ,  (32) 

11(4  -  ioio)zz"(4  -  Mo) ip 


where  Z  =  BZ  and  Bs  =  1„.  Let  R  be  the  true  covariance  matrix  associated  with  any  column  of  Z. 

Let  R  =  LL"  be  the  Cholesky  decomposition  of  R  where  L  is  an  A  x  LTM.  Thus  we  can  construct 

the  transformed  data  matrix  as  Z  =  LT  where  VisaN  x  K  matrix  of  i.i.d.  zero-mean  complex  circular 
Gaussian  random  variables  with  variance  equal  to  one.  Hence 


Now 


|zz"|  =  I|lvv“l"||  =  ||LiMI^^"l  ■ 


L,  Oh 

4^n  ' 

1 

o 

I 

_  0]^ 

(33) 


(34) 


where  L,  =  (A  -  1)  x  (A  -  1)  matrix  of  the  first  A  -  1  rows  and  columns  of  L,  y„  =  (A  -  1)  X  A 
matrix  of  the  first  A  -  1  rows  of  V,  =  last  row  of  V,  and  0„  denotes  a  column  vector  of  n  zeroes. 
Thus 


||(/,-M:)Zz"(/,-m:)||,=  l4^nV.‘/4"ll  =  WmMW 


(35) 


and 
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(36) 


/3  = 


II  I 

T^' 


If  4  =  (4«)>  II^P/  l|j^i  11^  =  Cw  Let  the  Q  -  R  decomposition  of  K  =  where  L,  is 

an  X  At  LTM  and  .  It  is  straightforward  to  show 


I  yy«  II 

^..Wi'll 


L^AWJ  J 


(37) 


where  Umf  is  the  sum  of  the  squared  magnitudes  of  AT  —  At  +  1  zero-mean  complex  circular  Gaussian 

r.v.’s  with  variance  equal  one  [11].  This  results  in  having  a  chi-square  pdf  of  order  2(K  -  N  + 

1).  Set 


/3  = 


"NN 


,M2 

•■NN  • 


Thus  jS  has  the  following  pdf 


1  1 

(K  -  N)\  p 

^  Ktl 


4 


K~N 

exp 


/3 


and 


P,  = 


(38) 


(39) 


(40) 


We  observe  that  the  MGLRT  for  a  range  distributed  target(s)  is  not  a  CFAR  detector  with  respect 
to  the  external  noise  environment;  i.e.  the  test  statistic  and  hence  the  false  alarm  probability,  Pp,  is  a 
function  of  a  parameter  of  the  true  covariance  matrix  of  external  noise  environment.  From  (39),  we  see 

that  the  test  statistic’s  pdf  is  a  function  only  of  4a/  which  depends  on  the  true  covariance  matrix  R  which 
in  turn  depends  on  the  true  covariance  matrix  R . 


Because  the  Pp  is  a  function  of  a  parameter  of  the  external  noise  environment,  there  would  seem  to 
be  no  way  of  a  priori  specifying  a  threshold,  T,  in  order  to  attain  a  given  Pp.  A  methodology  to  overcome 
this  deficiency  and  make  the  derived  adaptive  detection  scheme  practical  is  as  follows.  Let  the  first  NpQp 

columns  of  L  approximately  span  the  vector  space  associated  with  the  external  jamming  where  Np^p 
represents  the  number  of  degrees  of  freedom  (DOF)  necessary  to  span  this  space  and  Nppp  <  N.  The 
amplitude  levels  of  the  elements  of  the  first  Nppp  columns  tend  to  be  above  the  average  internal  noise 
amplitude  level  which  we  assume  without  loss  of  generality  equals  one.  The  last  N  —  Nppp  columns  of 

L  approximately  span  the  null  space  of  jamming  or  equivalently  the  internal  noise  space.  The  amplitude 
levels  of  elements  of  the  last  N  —  Nppp  columns  are  approximately  at  the  averge  internal  noise  amplitude 

level.  Hence  4a/  ~  L  the  jamming  scenarios  'when  Nppp  <  N,  it  will  be  found  that  the  external 
jamming  can  suppressed  at  or  below  the  internal  noise  by  the  non-adaptive  matched  filter  (i.e.  form 

s»R’‘Z). 
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A  methodology  for  insuring  that  a  given  Pp  is  not  exceeded  would  be  to  first  find  the  maximum  of 

over  all  the  expected  external  noise  scenarios  (in  a  systems  design,  this  is  normally  known).  If  the 
non-adaptive  matched  filter  is  working  well  (i.e.  suppressing  the  external  noise  below  internal  noise, 

which  in  most  cases  is  the  design  goal  and  drives  the  design  specification),  then  max  should  be  at 
or  about  the  internal  noise  level  (which  is  usually  normalized  to  one).  We  form  the  following  test  statistic: 


^  = 


/ 


NN 


max 


/(v)2  ^ 
‘AW 


(41) 


Thus 


Pp  =  Prob  {&  >  t\H^}  <  Prob  {/3  >  t\H,}  =  Pp.  (42) 

If  we  choose  t  based  on  obtaining  a  desirable  Pp,  then  the  actual  Pp  is  always  upper  bounded  by  Pp. 

Finally,  in  Section  IV  we  introduce  a  technique  which  is  CFAR  but  will  require  more  than  one  block 
of  input  data. 

B.  Probability  of  Detection 

In  most  cases  of  interest  the  probability  of  detection  will  be  computed  using  Monte  Carlo 
methods.  One  example  of  where  we  can  find  the  P^  explicitly  is  the  following.  Assume  that  range 
distributed  complex  amplitudes  of  the  desired  target  are  i.i.d.  zero-mean  complex  circular  Gaussian  r.v.’s 

with  variance  equal  to  a].  In  this  case  we  find  that  (see  Eq.  (37))  is  the  sum  of  the  squared 

magnitudes  oi  K  -  N  +  1  zero-mean  complex  circular  Gaussian  r.v.’s  with  variance  equal  to 

aj(s“P'‘s)  +  1,  [11].  Set  a],  =  aJ(s"P  ‘s).  We  find  that  for  the  test  statistic  given  by  (38)  that 

2((rjs  +  l)/w)J^  has  the  chi-square  pdf  of  order  2{K  -  N  +  \).  Hence  we  can  write  the  pdf  of  ^  under  the 
Hi  hypothesis  as 


P,(P\Hi) 


iliol  ^  1) 


K-N 

exp 

4(4  - 1) 

4(4  + 1) 


and 


P^  =  |4,(/3|H,)rf/3. 

IV.  A  CFAR  DETECTOR  FOR  RANGE  DISTRIBUTED  TARGETS 


(43) 


(44) 


From  the  test  statistic  form  given  by  (38),  it  is  seen  that  for  a  given  block  of  data,  Z,  that  the  test 
statistic  under  Hq  is  the  product  of  a  term  that  is  a  function  of  the  noise  environment,  4  that 

is  not,  .  We  can  create  a  CFAR  detector  by  using  the  MGLRT  test  statistics  over  multiple  blocks  of 
data  where  each  block  of  data  is  associated  with  non-overlapping  range  intervals.  We  assume  that  the 
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interferences  from  the  various  range  intervals  are  independent  but  that  any  snap  shot  of  data  (or  //-length 
vector)  share  the  same  N  x  N  covariance  matrix,  R. 

We  define  Z  to  be  the  primary  N  x  K  data  matrix  and  m  =  1,2,  . . .,  to  be  Mq  secondary 

N  X  K  data  matrices  where  we  assume  the  noise  vector  component  of  the  columns  of  Z  and 

ZJjn  =  1,  2,  ...,  Mq)  are  i.i.d  zero-mean  complex  circular  Gaussian  r.v.’s.  Let  ••• 

the  GLRT  test  statistics  for  the  primary  and  secondary  data,  respectively.  We  form  the  test  statistic 

^  K  (45) 

E 

/»  =  i 

If  we  assume  that  only  the  primary  data  can  contain  the  desired  signal,  then  (45)  takes  on  the  appearance 
of  a  cell-averaging  CFAR.  Based  on  the  form  of  the  individual  MGLRT  test  statistics  given  by  (28)  it 

is  seen  that  under  Hq  that  the  common  factor  of  cancels  from  the  numerator  and  denominator  of  (45) 
and  that  rj  is  independent  of  the  external  noise  environment.  Hence  (45)  defines  a  CFAR  detector  for  a 
range  distributed  target. 

We  point  out  that  we  could  have  also  formed  a  CFAR  detector  under  these  assumptions  using  the 
single  point  target  Kelly  detector.  For  this  we  would  perform  single  point  adaptive  target  detection  for 
each  range  cell  of  the  primary  data  and  use  the  secondary  data  (i.e.,  the  KMq  data  vectors)  as  the 
secondary  data  vectors  of  the  Kelly  detector.  After  doing  this,  a  1  out  of  K  detector  could  be  used  to 
determine  if  a  target  is  present.  A  disadvantage  of  using  this  detector  term  versus  the  proposed  CFAR 
test  (Eq.  (45))  is  the  computational  cost:  a  detection  statistic  must  be  computed  for  each  range  cell  using 
the  single  point  Kelly  detection  scheme  (this  includes  all  possible  ranges)  whereas  our  proposed  CFAR 
test  processes  the  ranges  in  blocks  of  K  ranges. 

Since  jS,  13 Jim  =  1,2,  . . .,  M^,)  are  all  chi-square  distributed  under  Hq,  it  is  clear  that  t]  will  have 
an  F  distribution  under  Hq.  In  fact,  if  we  set 


I3/[2(K  -  A  +  1)] 

15^  ’  (46) 

Y:  ^j[2Mq(k  -  a  + 1)1 

m=l 


then  pjn  I  Hq)  =  f(v;  2(K  -  N  *  1),  2Mq(K  -  M  +  1)),  where /is  the  pdf  associated  with  the  central 
F  distribution  with  parameters  p  =  2{K  -  A  +  1) ,  ^  =  2MJ,K  -  N  +  1) .  Thus 


Pf=  I  “  fir,  2{K  -  N  ^  1),  2M^{K  -  N  +  l))dr, .  (47) 

We  can  obtain  the  for  the  special  case  where  the  range  distributed  complex  amplitudes  of  the 
desired  target  are  i.i.d.  zero-mean  complex  circular  Gaussian  r.v.’s  with  variance  equal  to  a] .  It  will  be 
found  that 
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p^ir,  I  H,)  =  /  _J_;  2{K  -  N  ^  1),  2M„(/^  -  +  1)  ,  (48) 

<^«  +  1  I  +  1 

and 

Po  =  f  ”/  2(^:  -  +  1),  2M,(K  -  A^  +  1)  rfr,,  (49) 

I  +  1 

where  al  =  aj(s"^''s). 

V.  MGLRT  FOR  SIGNALS  ON  A  SUBSPACE 

In  this  section,  we  derive  the  MGLRT  for  a  range  distributed  target  or  targets  when  the  desired 
steering  vector,  s,  is  not  known  explicitly  but  is  defined  to  be  on  a  known  subspace  spanned  by  the 
columns  of  the  N  y.  matrix,  where  the  columns  of  are  orthonormal  and  2  ■<  ■<  N.  Thus 

=  1^  and  s  =  PjS  where  s  is  an  A4  x  1  vector  of  unknown  complex  values  with  the  constraint 
s”s  =  s“s  =  1. 

The  MGLRT  that  results  can  be  derived  in  similar  fashion  as  the  one  derived  in  Section  II  when  s 
is  assumed  known.  The  test  statistic  is  given  by 

sa'TlI^ 

with  constraint  s‘*s  =  1.  From  Section  II,  we  know  that 

Z“s  =  arg  min  ||(Z  -  sa'’)(Z  -  sa“)‘^||^.  (51) 

a 

Thus  we  desire  to  find 

§  =  arg  min  ||(4  -  P,ss«?")ZZ«(/^  -  P,ss«P,")»||^.  (52) 

s 

For  notational  convenience,  we  define  to  equal  the  product  of  the  positive  singular  values 

of  i4  and  recognize  that  ||A||^„  =  ||/1A”||^.  Thus  (52)  can  be  rewritten  as  the  equivalent  problem. 

s  =  arg  min  ||(4  -  P/s»P")Z||p,,-  (53) 

S 

Let  be  ihtNxN  LTM  associated  with  the  Q  -  decomposition  of  Z.  Then 

s  =  arg  min  1(7^  -  P,ss»P»)7:J|,„.  (54) 
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We  set  iFJ  where  is  an  N  x  (N  -  N)  matrix  whose  columns  are  the  orthonormal 

vectors  that  span  the  null  space  of  Thus  P^  is  an  x  unitary  matrix  and  P"P^  =  Now 


11(7^  -  P/s«Pr)i<JU  = 


.  SsHnH., 


^N-N, 

H 

Lr  - 

%-N, 

L 

sz 

s 

s 

(55) 


where  L„  is  the  LTM  associated  with  the  0  -  ^  decomposition  of  P^L  ■  Also 


- 

*■ 

r  “1 

\-N,  ' 

H 

0 

^SZ 

I.  - 

L  = 

N 

S 

S 

. 

SZ 

-  ss«)4f 

(56) 


where  we  have  written  as 


L 


(H) 

^N-N„N, 

(21) 

'sz 

(57) 


is  an  (N  -  NJ  x  (N  -  N,)  matrix,  '  is  a  A,  x  (77  -  N)  matrix,  and  is  a  A,  x  A,  matrix. 


It  is  seen  that  the  product  of  the  positive  singular  values  of  the  matrix  given  by  (56)  is  equal  to  the 
product  of  the  diagonal  elements  of  (these  elements  are  positive  with  probability  one)  and  the 
product  of  the  positive  singular  values  of  (/^  -  ss”)L^^’ .  Hence  the  minimization  of  the  far  RHS  of  (55) 
is  equivalent  to 


s  =  arg  mm 

S 


11(4  -  ss“)4 


(22) 

'sz 


psv 


(58) 


Write  the  SVD  of  as 


p(22)  _  p(22)  .(22)^(22)H 
J-'sz  ^sz  ^sz  ’ 


(59) 


where  and  are  A^  x  A^  unitary  matrices  and  is  an  A^  x  A^  diagonal  matrix  whose 
elements  are  the  A  singular  values  of  It  is  known  [1]  that  (58)  is  minimized  when  s  equals  the 
column  of  associated  with  the  largest  singular  value  of  Call  this  solution, 
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We  can  rewrite  the  MGLRT  test  statistic  as 


P  = 


(60) 


Let  a^,  n  =  I,  . . .,  be  the  singular  values  of  where  a^  >  >  ...  ^  <r^  andsetL^^  =  (1^^). 

Then 


n/<«’ 

^nun 

m=N-N^^\ 

li 


(61) 


The  procedure  for  finding  the  MGLRT  test  statistic  of  a  range  distributed  target  or  targets  with  s  on 
a  subspace  is  as  follows. 


1.  Set  =  [P,,  PJ 

2.  Find  the  N  x  N  LTM,  L^,  associated  with  Q  -  R  decomposition  of  Z. 

3.  Find  the  N  x  N  LTM,  associated  with  the  Q  -  R  decomposition  of  P^L^. 

4.  Find  the  singular  values  of  defined  by  (57). 

5.  Form  the  test  statistic  as  given  by  (61). 

We  point  out  that  the  distribution  of  /3  under  Hq  for  >  2  is  a  function  of  the  covariance  matrix  and 
2-thus  the  detector  is  not  CFAR.  Furthermore,  we  know  of  no  methodology  at  this  time  similar  for  the 
=  1  case  presented  in  Section  III  which  results  in  a  detector  which  is  CFAR  for  >  2. 

VI.  RESULTS 

In  this  section,  we  present  some  simulation  results  demonstrating  the  effectiveness  of  the  derived 
MGLRT  in  detecting  a  range  distributed  target  or  targets.  We  compare  the  performance  of  the  new 
detector  with  that  of  what  we  call  a  M  out  of  K  Kelly  detector  which  we  now  describe.  The  M  out  of  K 
Kelly  detector  consists  of  two  levels  of  detection.  The  first  level  is  a  single  point  scatterer  Kelly  detector 
[2]:  a  given  range  cell  in  theN  x  K  data  block  is  considered  the  candidate  or  primary  data.  This  primary 
data  is  an  A/^  X  1  vector.  The  other  K  —  I  range  cell  data  are  used  as  the  secondary  data  which  are 
assumed  in  the  derivation  of  the  Kelly  detector  not  to  contain  any  of  the  desired  signal.  This  secondary 
data  is  represented  as  a  A1  x  (K  —  1)  secondary  data  matrix.  The  standard  Kelly  detection' statistic  is 
computed  for  each  candidate  range  cell  (a  total  of  K)  where  the  candidate  or  primary  range  cell  is  varied 
over  all  possible  K  range  cells.  The  output  from  each  of  the  K  h)^othesis  tests  is  either  a  zero  or  one. 
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The  second  level  of  detection  consists  of  counting  the  number  of  detections  out  of  K  from  the  first  level 
detector  and  comparing  this  integer,  m,  to  the  integer  threshold,  M,  of  the  second  level  detector.  If  m  > 
M,  a  detection  of  a  target  in  the  N  x  K  data  block  is  declared.  A  disadvantage  of  the  M  out  K  Kelly 
detector  is  that  in  practice  we  must  set  M  =  1 .  Otherwise  for  scattering  scenarios  where  the  number  of 
scatterers  in  the  K  range  cells  is  less  than  M,  the  probability  of  detecting  the  target  on  targets  would  be 
almost  zero.  For  a  more  complete  description  and  analysis  of  the  M  out  of  K  Kelly  detector  see  [12]. 

We  present  results  for  a  single  representative  case  which  will  illustrate  the  performance  advantage 
that  the  MGLRT  of  a  distributed  target(s)  has  over  the  M  out  of  K  Kelly  detector  when  there  are  multiple 
scatterers  (or  targets)  in  range.  For  our  example  we  assume  the  inputs  to  the  adaptive  processor  are  the 
received  returns  for  an  N  element  linear  array  (half-wavelength  spacing)  with  no  time  taps.  We  use  the 
narrowband  assumption  so  that  the  desired  signal  vector  and  jamming  vectors  can  be  represented  as 

(1/\/A)  (1  .  Assume  a  single  jammer  is  the  only  source  of  external  interference  so 

that  the  W  X  covariance  matrix  is  given  by 

R  =  o-yJJ'^  +  (62) 

where  a]  is  the  power  of  the  jammer  referenced  to  the  internal  noise  power  of  each  antenna  element 

receiver,  J  =  (1  is  the  array  phasing  of  the  jammer,  and  4-  is  the  N 

X  N  identity  matrix  which  represents  the  covariance  matrix  of  the  internal  noise  contribution  of  each 
antenna  element  receiver. 

For  the  simulation,  we  choose  s  =  (1/\/A)  (1  1  ...  1)^,  the  number  of  samples  per  receive  channel, 
K  —  15,  the  number  of  sensor  inputs,  N  =  S,  and  the  probability  of  false  alarm  =  lO”'*.  In  Fig.  1,  the 
probability  of  detection  (Pp)  performance  of  a  single  point  scatterer  Kelly  detector  (not  the  M  out  K  Kelly 
detector)  and  the  MGLRT  of  a  distributed  target(s)  are  plotted  versus  the  input  signal-to-internal  noise 
power  ratio  for  the  case  when  there  is  only  one  scatterer  present  in  the  N  X  K  data  block  and  the 
jamming  power  is  zero.  It  is  seen  that  the  MGLRT  detector  of  distributed  targets  slightly  out-performs 
the  single  point  scatterer  Kelly  detector.  In  Fig.  2,  the  setup  is  the  same  as  Fig.  1  except  now  we  have 
two  distinct  (different  range  cells)  equi-powered  scatterers.  The  sum  of  the  SNR’s  of  each  scatterer  equals 
the  input  SNR.  We  see  clearly  that  the  single  point  scatterer  Kelly  detector  is  significantly  effected  by 
the  presence  of  a  desired  signal  (one  of  the  scatterers)  in  the  secondary  data.  In  fact,  as  the  power  of  the 
second  scatterer  inereases,  the  SNR  losses  become  much  greater  (this  was  also  shown  in  [7]).  We  show 
the  single  point  scatterer  Kelly  performance  in  order  to  demonstrate  that  the  M  out  of  K  Kelly  detector 
could  perform  poorly  in  a  multiple  scatterer  scenario  because  of  the  poor  detection  performance  in  the 
first  level  of  detection.  We  now  present  results  for  the  M  out  of  K  Kelly  detector  versus  the  GLRT  of 
distributed  target(s)  which  show  this. 

In  Fig.  3,  we  compare  the  Pp  performance  of  a  1  out  of  K  Kelly  detector  with  the  MGLRT  detector 
of  a  distributed  target(s)  of  a  single  scatterer  when  the  interference  scenario  consists  of  a  single  30  dB 
jammer  (with  angle  .6x)  and  internal  noise.  We  see  that  the  new  detector  slightly  out  performs  the  1  out 
of  K  Kelly  detector.  The  setup  is  the  same  for  Fig.  4  as  for  Fig.  3  except  two  distinct  equi-powered 
scatterers  are  present.  We  compare  the  1  out  of  K  and  2  out  of  K  Kelly  detectors  with  the  GLRT  of 
distributed  target(s).  The  2  out  of  K  Kelly  detector  is  superior  to  the  1  out  of  K  Kelly  detector  and  the 
new  detector  is  clearly  superior  to  the  2  out  of  K  Kelly  detector. 
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We  also  point  out  that  the  new  detector’s  detection  performance  is  about  the  same  for  when  there  are 
one  or  two  scatterers  present.  Also  the  detection  performance  is  about  the  same  for  the  example  when 
jamming  is  present  or  not.  For  the  jamming  scenario,  the  A^th  diagonal  element  of  the  Cholesky  factor 
of  R  is  approximately  equal  to  one.  Hence  based  on  the  results  of  Section  3A,  the  desired  thresholds  for 
a  given  for  the  two  noise  scenarios  are  approximately  equal. 

What  is  lost  by  using  the  MGLRT  of  a  distributed  target  or  target(s)  instead  of  a  single  point  scatterer 
Kelly  detector?  Clearly,  the  range  resolution  is  degraded.  With  the  new  detector  one  only  knows  that 
there  is  a  detection  among  the  K  range  cells.  If  we  are  looking  for  multiple  targets,  we  obviously  desire 
to  know  how  many  are  present  within  the  K  range  cells.  This  information  might  be  obtainable  if  after  a 
detector  indication  by  the  new  detector,  the  single  point  scatterer  Kelly  is  applied  to  each  range  cell  in 
order  to  see  if  individual  range  cells  are  detected  thus  restoring  to  some  degree  the  original  range 
resolution  capability  of  the  receiver. 

Finally,  the  new  detection  scheme  could  be  used  to  reduce  the  processing  load  of  a  range  cell-by 
range  cell  single  point  scatterer  adaptive  detector;  the  new  detector  would  be  used  as  a  first  level  detector 
indicating  whether  a  target  or  targets  are  present  in  a  set  of  range  cells.  A  second  level  of  detection  using 
the  single  point  scatterer  adaptive  detector  would  be  invoked  if  a  detection  recurred  in  the  first  level  of 
detection. 

VIL  A  COMMENT  ON  THE  MGLRT  FOR  A  SINGLE  POINT  TARGET 

Almost  all  of  the  adaptive  detection  literature  has  been  concerned  with  detection  of  a  single  point 
target;  i.e.,  the  target  is  modeled  as  a  single  point  scatterer  that  is  completely  contained  within  a  single 
range  cell.  The  data  that  can  contain  the  desired  target  returns  is  called  the  primary  data  vector.  Returns 
from  other  range  cells  are  used  to  form  secondary  data  vectors.  The  effects  of  contaminated  (by  desired 
signals)  secondary  data  on  performance  of  the  Kelly  detector  and  the  AMF  detector  were  examined  in 
[7]  for  the  case  where  the  target  is  modeled  as  a  point  scatterer.  It  was  shown  that  significant  degradation 
in  detection  performance  can  occur  when  multiple  targets  of  equal  scattering  strength  are  contained  within 
the  same  block  of  primary  and  secondary  data.  The  two  detectors,  the  Kelly  and  AMF,  do  not  take  into 
account  the  possibility  of  contamination  of  the  secondary  data  by  other  target  returns.  The  search  for  an 
effective  adaptive  detector  against  secondary  data  that  has  been  contaminated  has  been  a  long  standing 
unsolved  problem.  In  this  section  we  derive  the  MGLRT  for  this  problem.  We  assume  that  the  steering 
vectors  associated  with  the  contamination  are  identical.  In  practice,  this  is  often  the  case,  since  in  a  given 
range  interval  which  is  used  to  construct  the  primary  and  secondary  data,  multiple  targets  are  usually 
transiting  at  the  same  speed  and  direction  (such  as  a  formation  of  air  targets). 

Interestingly,  it  will  be  found  that  the  derived  MGLRT  does  not  use  the  secondary  data  at  all.  Hence 
the  derived  GLRT  is  not  useful  as  an  adaptive  detector  since  one  normally  uses  the  secondary  data  to 
estimate  the  unknown  covariance  matrix  of  the  received  interference.  As  a  result,  the  MGLRT  formalism 
for  solving  this  problem  is  not  the  solution  and  the  search  for  an  effective  adaptive  detector  against 
secondary  data  is  still  an  open  problem. 

We  now  briefly  outline  the  derivation  of  the  MGLRT  for  single  point  target  and  show  that  the  test 
statistic  is  independent  of  the  secondary  data  set.  We  collect  data  from  sensors  and  form  a  primary  data 
vector  associated  with  a  given  range.  We  also  collect  data  from  the  N  sensors  at  AT  -  1  other  range  cells 
and  form  K  —  I  secondary  data  vectors  where  other  targets  returns  may  be  present.  The  hypotheses  are 
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H,: 


^0  =  ^0 


=  x^.  +  a/s  k  =  1,  2,  K  -  1 


(63) 


H,: 


=  x^,  +  a/s  k  =  I,  2,  ...,  K  -  I , 


where  Zq  is  the  primary  data,  z^(k  =  I,  2,  K)  are  the  secondary  data  vectors,  Xq  and 

x^(k  =  I,  2,  K)  are  the  interferences  and  a^ik^O,  1,  K  -  1)  are  the  complex  signal 
amplitude  associated  with  the  A:’th  range  cell.  Similarly  to  the  development  of  Section  II,  (with  the  same 
assumptions  on  the  interference),  we  can  derive  the  MGLRT  statistic  as 


min  II (Z^  -  sao)(Z^  -  +  Vo”  11/ 

“o _ 

min  |1(Z  -  sq:“)(Z  -  sa/VIL 


(64) 


where  Z^  =  (z,  z^ . . .  z^_,)  is  the  N  x  (K  -  1)  matrix  of  secondary  data,  Z  =  (z^  i  Z)  is  the  Nx  K  total 
data  matrix,  =  (a,,  a^,  ...,  and  a,  =  (a,,  a,,  ...,  a^_,y.  Let  be  a  unitary  matrix  such 

that  QjZo  =  1^0 1^0  where  |zj  =  Define  Z^  =  Q^Z^,  Z  =  Q^Z  =  (Izjlo  i  Z),  and 

s„  =  (2  s.  Thus  (64)  can  be  rewritten  as 

0  Zq 


/3  = 


min  II (Z^  -  s^ao)(Z^  -  s^otoV^  +  |Zol%lo^l 

"o  _ 

min  II  (Z  -  Sga:")(Z  -  Soa,“)||^ 


(65) 


The  solution  for  a,  was  derived  in  Section  II: 


a,  =  SoZ. 


(66) 


The  solution  for  ao  in  the  numerator  of  (65)  is  now  considered.  The  numerator  can  be  rewritten  as 

NUM  =  min  ||(Z  -  SoOf“)(Z  -  SQa")“||^,  (67) 


where  a  =  (0  ;  a^)^.  It  is  straightforward  to  show  that  there  exists  an  a  of  this  form  such  that  Z  - 
has  rank  N  -  1.  Hence,  we  look  for  solutions  for  a  on  the  space  where  Z  -  s^a”  has  rank  -  1. 

We  will  need  the  following  lemma  in  order  to  proceed. 

Lemma  2:  If  Z  -  SgO:”  has  rank  N  -  1  and  Z  has  rank  N,  then  there  exists  an  JV-Iength  vector  s,,  such 
that 


<Jn  -  SoSo)Z  =  Z  -  s^a". 


(68) 
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with  s"sg  =  1. 

proof:  See  Appendix  B  □ 

Using  Lemma  2,  we  see  that  finding  a  is  equivalent  to  finding 

§0  =  arg  min  ||(/^  -  s„So”)Z||^^^.  (69) 

So 

Set  Sq  =  (Sq,  §02  §03  ...  §0^)^.  Because  a”  =  §“Z  =  §^\|Zo|  lo  ;  Z)  and  the  first  element  of  a  is  zero, 
it  follows  that  Jq”,  I^oI  =0.  Since  |Zo|  5^  0  with  probability  one,  the  first  element  of  §0  must  be  zero 
or  ioi  =  0  which  is  a  constraint  on  the  minimization  problem  posed  by  (69)  along  with  the  constraint 

So\  =  1.  Using  a  methodology  similar  to  that  of  finding  the  solution  for  (23),  we  can  show  that  (69) 
is  equivalent  to  solving  the  following  problem 

min  §"§Q  ,  constraints  s^,  =  0  ,  §"Sq  =  1 .  (70) 

s 

It  is  straightforward  to  show  that  the  optimal  solution  is  given  by 


(71) 


(72) 


(73) 


The  elements  of  are  derived  from  the  elements  of  Sq  which  equals  Q^.  Thus  Sq  is  only  a  function  of 

the  primary  data  Zq  and  not  a  function  of  any  of  the  secondary  data,  Zj.  Hence  this  detector  is  not  useful 
since  it  essentially  throws  away  the  data  that  would  be  used  to  estimate  R. 


Recall  for  the  adaptive  detector  of  a  range  distributed  target  that  under  7/,  (signal  plus  noise)  that  the 
MGLRT  projected  the  data  matrix,  Z,  onto  the  null  space  of  the  desired  signal  vector.  Hence,  a  hueristic 
adaptive  detection  scheme  for  a  point  target  with  contaminated  secondary  data  would  be  to  project  the 
secondary  data  matrix  onto  the  null  space  of  s  under  the  hypothesis.  The  resultant  heuristic  detector 
has  the  form 
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(74) 


_  i(/^  -  -  ss“)  +  zX 

11(4  -  ss")ZZ"(4  - 

This  also  turns  out  not  to  be  an  effective  detector.  In  fact,  the  detection  statistic  is  approximately  equal 
to  |s”Zg|^  where  again  the  secondary  data  is  not  used. 

VIII.  CONCLUSIONS 

A  Modified  Generalized  Likelihood  Ratio  Test  (MGLRT)  for  the  adaptive  detection  of  a  target  or 
targets  that  are  distributed  in  range  was  derived.  The  unknown  parameters  associated  with  the  hypothesis 
test  are  the  complex  amplitudes  in  range  of  desired  target  and  the  unknown  covariance  matrix  of  the 
additive  interference  which  is  assumed  to  be  characterized  as  complex  zero-mean  correlated  Gaussian 
random  variables.  The  target’s  or  targets’  complex  amplitudes  are  assumed  to  be  distributed  across  the 
entire  input  data  block  (sensor  x  range).  Results  on  probabilities  of  false  alarm  and  detection  were 
derived.  The  MGLRT  statistic  was  found  to  be  dependent  on  the  external  noise  environment  and  thus 
does  not  have  the  desirable  quality  of  being  a  CFAR  detector. 

Methods  were  presented  which  would  make  it  almost  CFAR  (upper-bounded  if  only  one  input 
data  block  were  used  and  CFAR  if  multiple  input  data  blocks  were  used.  Some  simulation  results  were 
presented.  It  was  shown  that  the  derived  MGLRT  of  range  distributed  targets  is  much  more  effective  in 
detecting  targets  distributed  in  range  than  a  M  out  of  K  detector  cascaded  with  a  single  point  target  Kelly 
detector.  In  fact  a  solution  to  abate  the  deleterious  effect  of  the  secondary  data  being  contaminated  by 
desired  signal  is  to  detect  all  the  target  scatterers  within  the  data  block  using  the  MGLRT  for  range 
distributed  targets  (such  as  a  formation  of  targets). 

Finally,  the  MGLRT  associated  with  detecting  a  single  point  target  with  signal  contaminated 
secondary  data  was  also  derived.  It  was  shown  that  surprisingly,  this  detector  is  ineffective  in  solving  this 
problem. 
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Appendix  A 

PROOF  OF  LEMMA  A1 


Lemma  Al:  Let  A  h&  a.  K  x  N  matrix,  and  C/,,  U2  be  two  N  x  M  matrices  such  that 

1.  Sp{U\)  =  Sp{U^  where  Sp  denotes  Span  of, 

2.  M  <  N, 

3.  U^U.  =  /^,  /  =  1,  2, 

then 

proof:  Since  Sp{U\)  =  Sp{U^,  we  can  find  anM  x  M  matrix  B  such  that  f/,  =  Now 

=  /^.  (Al) 

Thus  B  is  unitary  and 

'  UUjiW^  =  lum,,..  (A2) 

□ 
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Appendix  B 

PROOF  OF  LEMMA  2 


Before  proofing  Lemma  2,  we  prove  the  following. 

Lemma  BJ:  Let  Ahea.  N  X  N  nonsingular  matrix,  s  an  A  x  1  vector,  and  b  a  AT  x  1  vector  with  K  > 
N.  Define  B  =  [A  :  -  sb".  If  B  has  rank  N  -  I,  then  the  N  +  1st  through  ATth  elements  of  b 

are  zero. 

proof:  The  matrix  BB"  also  has  rank  N  -  1.  We  can  write 

=  (A  -  sbo“)(A  -  sb")  +  ass",  (Bl) 

where  bo  is  a  A-length  vector  equal  to  the  first  N  elements  of  b  and  a  equals  the  sum  of  the  magnitude 
squares  of  the  last  K  -  N  elements  of  b.  Now  if  B  has  rank  N  -  1  there  exists  an  A^-length  vector,  u, 
such  that  u"5  =  0.  Thus 

=  u"(/4  -  sb")(^  -  sbo”)"u  +  a|u"s|^  =  0.  (B2) 

Since  {A  -  sbo")(/4  -  sbo")"  is  a  positive  semi-definite  matrix,  (B2)  implies  if  a  5^  0  that  u"s  =  0.  In 
this  case  it  follows  that  u"AA"u  =  0  which  contradicts  the  assumption  that  A  is  nonsingular.  Therefore 
a  =  0  and  Lemma  B1  follows  ^ 

We  now  prove  Lemma  2.  Let  Z  =  [L 0^^k-n\Q  be  the  Q  -  R  decomposition  of  Z  where  L  is  an  A 
X  N  LTM  and  Q  is  a  a:  X  a:  unitary  matrix.  Now  if  Z  -  sa"  has  rank  N  -  1  then  [L :  0^  J  -  sa^ 
has  rank  N  -  1  where  •  From  Lemma  Bl,  the  last  K  -  A  elements  of  ao  must  be  zero.  Thus 

there  exists  a  s  such  that 


ao"  =s»L, 

where  ao^  are  the  first  N  elements  of  ao  and  L  is  nonsingular.  Hence 

a"  =  s“[L:0]e  =  s«Z, 


and  the  lemma  follows. 


(B3) 


(B4) 

□ 
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Probability  of  Detection  (Pd) 


Two  Scatterers,  Kelly  (solid),  New  Detector  (dashed),  N=8,  K=15,  Pf=.0001 ,  white  noise 
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Fig.  1  —  Performance  comparison  of  a  single  point  scatterer  Kelly  detector  with  the  GLRT  detector  of 
a  distributed  target(s):  one  scatterer,  N  =  8,  K  =  15,  =  10"^  white  noise  interference. 


2. 


Probability  of  Detection  (Pd) 


Single  Scatterer,  Kelly  (solid),  New  Detector  (dashed),  N=8,  K=15,  Pf=.0001,  white  noise 


Fig.  2  —  Performance  comparison  of  a  single  point  scatterer  Kelly  detector  with  the  GLRT  detector  of 
a  distributed  target(s):  two  scatterer,  N  =  S,  K  =  15,  =  10■^  white  noise  interference. 
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Probability  of  Detection  (Pd) 


Probability  of  Detection  (Pd) 


Two  Scatterers,  N=8,  K=15,  Pf=.0001 ,  J/N=30  dB 


Input  SNR  (dB) 

Fig.  4  —  Performance  comparison  of  a  1  out  of  K  Kelly  detector,  2  out  of  K  Kelly  detector  with  the 
GLRT  detector  of  a  distributed  target(s);  two  scatterers,  N  =  ^,  K  =  15,  Pp  =  10"^,  interference  is  a 
30  dB  jammer  plus  white  noise. 
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